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We consider small-angle electron-positron scattering in Quantum Electrodynamics. 
Leading logarithmic contributions to the cross-section are explicitly calculated to three 

* i 

loop. Next -to-leading terms are exactly computed to two loop. All the radiative cor- 
rections due to photons as well as pair production are taken into account. The impact 
of newly evaluated next-to-leading and higher order leading corrections is discussed and 
numerical results are explicitly given. The results obtained are generally valid for high 
and low energy e + e~ colliders. At LEP and SLC these results can be used to reduce the 
uncertainty on the cross-section below the per mille level. 
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An high accuracy, better than one per mille, measurement of the luminosity has been 
reached at LEP The small-angle electron-positron scattering (Bhabha) process is normally 
used as the reference cross-section to measure the luminosity. 

A poorly known Bhabha cross-section has the consequence of producing a systematic error 
on the determination of relevant physical observables as, for example, the hadronic peak cross- 
section (Jp eak or the leptonic widths T e ^. This poor accuracy reflects itself on the extraction 
of the Standard Model parameters on the exclusion of new physics signals as, for example, 
those coming from a deviation from the Standard Model prediction for the number of light 
neutrinos 0. The experimental uncertainty in the luminosity determination together with 
the theoretical one define the uncertainties on the above mentioned quantities. An adequate 
theoretical accuracy on the cross-section is therefore highly needed. 

Various approaches have been used to obtain the Bhabha cross-section || |J. Some of them 
(see ref. ||) are based on Monte Carlo generators. Others use the structure function method 
to calculate the radiative corrections ||. With both methods the contributions of next-to- 
leading as well as higher order radiative corrections have not been systematically evaluated 
and this yielding an important source of uncertainty above the level of 8a jo ~ 0.001 accuracy. 
Therefore an equally accurate theoretical determination of the Bhabha cross-section, even if 
it is approaching to, has not yet reached the experimental accuracy @. 

In this letter we describe the results obtained with a different approach based on the 
direct evaluation of Feynman diagrams. The results, given in analytical form, systematically 
take into account leading as well as next-to-leading contributions thus reducing the physical 
uncertainty to the 0.1% level. 

The differential cross-section for the small-angle Bhabha process in the Standard Model 
is §: 

= ^(l + * + *weak). 6 « 1 ' 8' = -° 2 / 2 + 9 ^/40, (1) 

5 weak = 2 9v^-^(9v+9 2 aWex + —(gt + 9t + Q9v9 2 a)\x\ 2 , 

X = s[(s-M 2 z + iM z T z )sm(29 w )}-\ £ = t[(t - Mf) sin^)]" 1 , 

where 9 is the scattering angle, s = 4E 2 , t = —E 2 6 2 , E is the beam (center-of-mass) energy, 9w 
the Weinberg angle and g a = —1/2, g v = —(1 — 4 sin 2 0^0/2 the axial and vector couplings of 
the Z° boson. Experimental cross-sections are obtained by collecting events within particular 
angular portions of the detectors with given energy cuts. To compare with those observed 
distributions one has to take into account radiative corrections to higher orders in perturbation 
theory and it is convenient to introduce the dimensionless quantity: 

O 2 ) ) r r A a e+e-^ + {q^,x 2 )e-{qi,xi)+x 
" 1 ' dxi / dx 2 e( Xl x 2 - x c ) / dV@i / dV©2 ; — in i 2 i , (2) 



" " 47^ J ^ J — ^ J - « -1 J - « ^ dxid 2 q j- dx2d 2 q ± 

where x c is the energy fraction threshold for the detection of the final electron and positron 
x c < x±x 2 , where xx >2 are the energy fractions of the scattered leptons and q^, q^r are the 
components of their momenta transverse with respect to the beam directions. 0^ 2 are the step 
functions which account for the registration of the scattered leptons by the circular detectors 
situated close to the beams with an aperture bounded by the angles #i,2,3,4 : 

e.-efc-MjraML*), e5 = e ( fc-Jg)e ( Jg-fc). (3) 
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If 6\ is the minimal aperture angle we define the minimal momentum transferred from the 
electron to the positron as Q\ = E 2 6\. At LEP/SLC Q\ ~ 1 GeV/c. For symmetric detectors: 
9i = 6 2 , 63 = 64, p = 63/61. 

E = S + E 7 + E 77 + E 7 + E e+e ~ + E 37 + E e+e ~ 7 (4) 
= E°(l + 8 tot ) = E°(l + 5 + 5 7 + 5 27 + 5 e+e " + 5 37 + cf +e ~ 7 ) 
= E°(l + ^5 i ), 5^ = 5^ + 5^, E° = l-i 



/' 2 



E has the form: 



dz 



So = / ^tf(-*g?)(i + <W(* 2 ) + W 2 )) 



(5) 



with £/(£) = (1 - n(t))" 2 where II(t) is the vacuum polarization operator of the exchanged 
photon ||. The experimentally observable cross-section is given by the expression: 



Ana' 



a 



Q 



E°(l + S + <5 7 + 5 27 + 6 e+e ~ + 5 37 + <5 e+e ~ 7 ). (6) 



The quantity E 7 collects the radiative corrections to the Born amplitude related to the 
emission of single virtual, real soft or hard photon. Virtual and soft photon contributions to 
the differential cross-section are proportional to the Born amplitude: dcr = der^l + -(2(L — 
1) In A + 3L/2 — 2)). This quantity contains the large logarithm L, L = ln(zQ 2 /m 2 ) (m is the 
electron mass), and A = 5E/E, where 8E is the maximal energy carried by a soft photon, 
A <C 1. The A-dependence disappears in the total sum when the emission of a hard photon 
with energy fraction larger than A is also taken into account. 

A simple analytical expression for the small-angle Bhabha differential cross-section for 
one hard photon emission can be written within the the infinite-momentum frame formalism. 
One obtains: 

da H 2a 2 1 + x 2 \(q 2 L ) 2 (l-x) 2 2m 2 x(l - x) 2 {d 1 - d 2 ) 2 ] ^ 



d x d 2 {l+x 2 )d 2 d 2 



dxd 2 qtd 2 q 2 L ~ 7i 2 ((^) 2 ) 2 (1 - x 
d l = m 2 (l - xf + {q^ - q^f = (1 - x)2 Pl k, 
d 2 = m 2 (l - x) 2 + (q^ - xq^) 2 = x(l - x)2q 1 k, q^ = q± + k 1 , 

x = -j, , qi=xp 1 + q 1 , pi + p 2 = q% + k + q 2 , 
Pi 

where pi, p 2 and k are the 4- momenta of the initial electron, positron and emitted photon 
respectively, x and 1 — x are the energy fractions of the scattered electron, is the transverse 
component of the photon momentum with respect to the beam direction. The photon is 
supposed to be emitted from the electron line. The same result is obtained for the emission 
from the positron line. The above expression takes into account only scattering-type diagrams. 
According to the definition of eq.(^) we have that 

E 7 = - J J dxU(-zQj) UL - l)P(x)[l + Q(x 2 p 2 - z)] (8) 
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+ 



1 + x z 



k(x, z) 



x 

- X) 



k(x, z) — 5(1 — x)\, 



1 + x 2 



[1 + 6(xV - z)\ + L a + <3>(x 2 p 2 -z)L 2 + 0(z 



x 2 p 2 )L 3 



L x = In 

and 



x 2 (z-l)(p 2 -z) 



(x — z)(xp 2 — z) 



In 



(z — x 2 )(x 2 p 2 — z) 



x 2 (x — z)(xp 2 — z) 



In 



(z — x 2 )(xp 2 — z) 



(x — z)(x 2 p 2 — z) 



P(x) 



1 + x 2 
1 — X 



Hm {l^e(l-x-A) + (^ + 21nA) 6{l-x)} 



(9) 



is the kernel of the non-singlet (NS) evolution equation ||. 

For small scattering angles (— t/s ~ 9 2 <C 1) the elastic amplitude A(s, t) has a generalized 
eikonal form: A(s,t) = Ag orn (t) T 2 (t) Ju(t) exp{i<p(t)}, where T(t) is the Dirac electron 
form factor and (f>(t) is the Coulomb phase ||. The Coulomb phase factor collects the multiple 
exchange photon contributions in the scattering channel. This holds for elastic as well as 
for inelastic amplitudes. This form permits us to neglect diagrams with several exchanged 
photons since the phase factor exp{i<f)(t)} disappears in the physical cross-sections. To higher 
orders in the perturbative expansion, only Feynman diagrams with one exchanged photon 
between electron and positron line do contribute. For the second order photonic contributions 
we use the known results of the electron Dirac form factor [10). Pauli form factor gives a 
negligible contributions of 0(8 2 (a/iT) 2 ). They include the cross-section for the emission of 
two soft photons as well as the emission a single real (soft or hard) photon with one-loop 



virtual radiative corrections |Tl| and the cross-section for the emission of two hard photons 
in the same or in opposite directions. For this last kinematical configuration we may use a 
generalization of the single bremsstrahlung cross-section in eq. ( 
expression: 



12| . We obtain the following 



^7 



4 7T 



J^U(-Q\z){l 2 J d Xl J dx 2 P(x 1 )P(x 2 ) 



(10) 



X C /X1 



x [Q{z - l)Q{p 2 -z) + Q{z - x\)Q{xlp 2 - z)\ 
x [&(z - l)0(p 2 -z) + &(z - x 2 2 )Q(x 2 2 p 2 - z)\ + L&X 



with $^ an analytical function (see |7]]). To calculate within the single logarithmic accuracy 
the contribution of the emission of two hard photons from a single lepton line (the electon 
in our case) we separate, by introducing a new auxiliary parameter 9 (1 ^> # ^> m/e), 
in collinear and semi-collinear kinematical regions ||13|| . Namely, in the collinear region both 
emitted photons move in the narrow cones defined by the polar angle 8q with respect to the 
initial or final (as well as simultaneously) electron 3-momenta. This region gives leading and 
next-to-leading contributions. The semi-collinear region produces only next-to-leading terms 
and corresponds to the kinematical configurations when only one of the photons moves inside 
one of the defined narrow cones when the second is radiated outside. It can be shown that, 
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as it should, the dependence on # disappears in the total sum: 



I(^ 2 

2 71 ' 



p(2)( 



X 



j^U(-Qlz)U 2 j dx {\p^(x) [ Q(x 2 p 2 -z) + l) 

1 ^ x c 

+ J jP(t) P(f) e(tV - z)} + L*rr\ 

X ' 



:n) 



lim 



21nA + 



+ 2 



dt 

n 

1+lVnl /, \ , 3 

2 ln(l - x) - In x + - 

l-x\ K ' 2 



5(1 -x) 



'1 + x) lnx — 1 + x 



(12) 



0(1 



re - A)}. 



For $ 77 see [0]. To the same order of perturbation theory we have to take into account the 
e + e _ pair production processes. For hard pair production we again consider four collinear 
kinematical regions, when the created pair moves close to the directions of the projectiles or 
of the scattered particles, and six semi-collinear regions which having 2^3 like kinematics 
i.e. when one of the final particles moves close to the beam direction or when one component 
of the created pair is close the other. The cancellation of the auxiliary parameter 8q can be 
explicitly verified ||14j| . 

The relevant two-loop contribution from the electron form factor [|l(J contains L 3 terms 
which disappear in the sum with the contribution due to soft pair production [JTJ]: 

i 



R{x) 



la, 2 

4 l 7T 

2 



dz . 



J -2 U (-zQi){ L J dx[l + Q(z - l)Q(x 2 p 2 -z)} R(x) + L$ e 

1 X c 

1 - x , . _ . 



(13) 



Fix) + 2(1 + x) In x + - (4 + 7x + Ax 2 

3 3x 



where $ e e collects nonleading corrections which can be found in J/J. 

Within the third order of perturbation theory it is sufficient, to the accuracy of the O(10~ 3 ), 
to consider only leading contributions related to the initial particle radiation. Leading contri- 
butions can be systematically included by using QED evolution equations ||. 

We consider the channels of 777 and 7e + e~ production real or virtual in all possible 
combinations. We obtain 

1 1 



S 3 7 



1 ,a 
4 n 



dz 



dxi 



dx-, 



Xc/Xl 



6 



5(1 - X2 ) p (3) (xi; 



(14) 



1 



x Q(x\p 2 - z)Q(z - 1) + ^P {2) (x l )P(x 2 )Q 1 Q 2 



e lj2 

P^(x) 



Q(z-x 12 )Q(p x l2 -z), 



y y 



£ = ln 



Ql 



1 ,a 
4V 



£) 1 



dz 



Xc/Xl 



d Xl I dx 2 \^[R p (x 1 )-^R s (x 1 )} 



(15) 
(16) 



x 5(1 - x 2 )Q(x{p 2 - z)Q(z - 1) + - P(x 2 )R r (x 1 ) 
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where 



R r (x) = R s (x) + -P(x), R s (x) = - — -U + 7x + Ax 2 ) + 2(1 + x) lnx, (17) 

3 3a; 



R r {x) = R s (x)\-+2\ii(l-x)J + (1 + x)[-\n A x - A J dy 

o ^ 

+ i(-9-3a; + 82; 2 )lna; + -f---8 + 8x + 3x 2 ) + -P {2) (x). 
3 3\ i / 3 

By combining the partial results in eqs. (10-14), and ([16]) one obtains the final result 
for the observable cross-section eq.(|6]). 

In Table 1 we give our results for different values of the threshold energy fraction x c for 
the defined angular acceptance of Q\ = 1.6° and 9 2 = 2.8°. Nonleading second order photonic 
corrections are negative and turn out to be larger in magnitude than both third order and 
second order ones due to pair production. This is shown in Fig. 1. Photonic corrections 
(leading and nonleading) (see Fig. 2) due to double photon emission from both fermions 
dominate the ones from a single lepton line. This fact is a consequence of the sign-changing 
p( 2 )(x) function (see eq.(|T^)), which describes double bremsstrahlung from a single lepton 
line. The importance of radiative corrections to Bhabha cross-section is illustrated in Fig. 3. 
The decreasing of the cross-section in the region x c — > 1 can be understood as a reduction 
by increasing x c of the positive contribution of real photon emission, while first order virtual 
corrections, being negative, remain unchanged. 

The accuracy of this result is implicitely defined by the terms omitted in the perturbative 
series. Typically they are of the type: 

-Q\ ft 2 , A 3 L 2 , (— L) 4 . (18) 

7T 7T 7T 7T 

An estimate of their magnitudes permits us to state that the uncertainty of our result to 
be at the level 10~ 4 . A detailed analysis of the attained accuracy is given in 0. From the 
analysis given above clearly emerges the importance of the next-to-leading contributions to 
reach an accuracy adequate to the one attained in present LEP and SLC experiments. The 
above formulae can be applied to future high-energy e + e~ colliders as well. 
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1 Figure Captions 



Fig.l: The ratios a) Sll nleading /5 to t, b) 5 e+e / 5 tot and c) (<5 37 + 5 e+e 7 )/5 tot as functions of x c 
for Q x = 1.6°, 6 2 = 2.8° and ^fs = 91.161. 

Fig.2: Behaviour of a) SH ding (solid line) and 5% ading + SH nleading (dashed line), b) &! lleading 
(solid line) and 52, .. + <G , ,. (dashed line) as functions of x c (parameters are as in 

v ' Heading 1 nonleading v ' vi 

Fig.l). 

Fig.3: Corrected cross-section of Bhabha scattering according to eq.(^) as a function of x c : a) 
<Jo+0(a) corrections, b) with 0(a 2 L 2 ) photonic corrections added, c) with all other corrections 
(see eq.(|)) added. a is the Born cross-section. For the parameters of as in Fig.l one has 
cr = 106.33 nb. 



Table 1: Per cent values of 5 l as defined in eq.(f|) for y/s = 91.161 GeV, Q\ = 1.61°, 
6 2 = 2.8°, sin 2 6 W = 0.2283, Y z = 2.4857 GeV. Quantity 5 2 ~< is equal to the sum 5^ + 5^ (see 
eq.(|)). 



x c 


So 




^leading 


u nonleading 


5 e+e ~ 




5 3 7 




0.1 


4.120 


-8.918 


0.657 


0.162 


-0.016 


-0.017 


-0.019 


-4.031±0.006 


0.2 


4.120 


-9.226 


0.636 


0.156 


-0.027 


-0.011 


-0.016 


-4.368±0.006 


0.3 


4.120 


-9.626 


0.615 


0.148 


-0.033 


-0.008 


-0.013 


-4.797±0.006 


0.4 


4.120 


-10.147 


0.586 


0.139 


-0.039 


-0.005 


-0.010 


-5.356±0.006 


0.5 


4.120 


-10.850 


0.539 


0.129 


-0.044 


-0.003 


-0.006 


-6.115±0.006 


0.6 


4.120 


-11.866 


0.437 


0.132 


-0.049 


-0.002 


-0.001 


-7.229±0.006 


0.7 


4.120 


-13.770 


0.379 


0.130 


-0.057 


-0.001 


0.005 


-9.194±0.006 


0.8 


4.120 


-17.423 


0.608 


0.089 


-0.069 


0.001 


0.013 


-12.661±0.006 


0.9 


4.120 


-25.269 


1.952 


-0.085 


-0.085 


0.005 


0.017 


-19.379±0.006 
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Figure 2: Behaviour of a) 5% ading (solid line) and 5% ading + ^nieading (dashed line), b) V' leadmg 
(solid line) and 52, „. + 52 , ,. (dashed line) as functions of x c (parameters are as in 

v ' 1 leading 1 nonleading v ' vi 

Fig.l). 
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Figure 3: Corrected cross-section of Bhabha scattering according to eq.(6) as a function of 
x c : a) o"o + @{ a ) corrections, b) with 0(a 2 L 2 ) photonic corrections added, c) with all other 
corrections added, where <r is the Born cross-section. For the parameters of as in Fig.l one 
has (T = 106.33 nb. 
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